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A numerical approach is considered for spherically symmetric spacetimes that generalize LemaTtre- 
Tolman-Bondi dust solutions to nonzero pressure ("LTB spacetimes"). We introduce quasi-local 
(QL) variables that are covariant LTB objects satisfying evolution equations of Friedman-Lemaitre- 
Robertson-Walker (FLRW) cosmologies. We prove rigorously that relative deviations of the local 
covariant scalars from the QL scalars are non-linear, gauge invariant and covariant perturbations 
on a FLRW formal "background" given by the QL scalars. The dynamics of LTB spacetimes is 
completely determined by the QL scalars and these exact perturbations. Since LTB spacetimes 
are compatible with a wide variety of "equations of state", either single fluids or mixtures, a large 
number of known solutions with dark matter and dark energy sources in a FLRW framework (or 
with linear perturbations) can be readily examined under idealized but non-trivial inhomogeneous 
conditions. Coordinate choices and initial conditions are derived for a numerical treatment of the 
perturbation equations, allowing us to study non-linear effects in a variety of phenomena, such as 
gravitational collapse, non-local effects, void formation, dark matter and dark energy couplings and 
particle creation. In particular, the embedding of inhomogeneous regions can be performed by a 
smooth matching with a suitable FLRW solution, thus generalizing the Newtonian "top hat" models 
that are widely used in astrophysical literature. As examples of the application of the formalism, we 
examine numerically the formation of a black hole in an expanding Chaplygin gas FLRW universe, 
as well as the evolution of density clumps and voids in an interactive mixture of cold dark matter 
and dark energy. 
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I. INTRODUCTION. 

Recent observations suggest that present cosmic dy- 
namics is dominated by elusive sources called "dark mat- 
ter" and "dark energy" , the former clustered in galac- 
tic halos and the latter possibly associated with a large 
scale repulsive (yet unknown) interaction. A large body 
of phenomenological and theoretical models have been 
proposed to describe these sources (see [T] for a compre- 
hensive review). The dominant model of dark matter is 
a collision-less gas of supersymmetric particles: "cold" 
dark matter pQ, whereas dark energy has been described 
by a cosmological constant [5] , as well as by quintessence 
scalar fields, Phantom fields, tachyons, branes (see [T]). A 
unified perspective is given by the Chaplygin gas [TJ|3j|4j, 
a single source that exhibits the expected dynamical be- 
havior of dark matter and dark energy. There is also a 
large number of empiric and phenomenological descrip- 
tions of dark matter interacting with dark energy [5j [6] 
(see [7] for a recent appraisal), as well as empiric dark 
energy "equations of state" that fit observations [8J|9]. 

While dark matter at the galactic scale is considered 
inhomogeneous, due to its association with structure for- 
mation, most research work on dark energy sources has 
been conducted in homogeneous Friedman-Lemaitre- 
Robertson- Walker (FLRW) spacetimes and/or linear 
perturbations, since dark energy is assumed to have a 
significant effect only at a larger cosmic scale in which 
the universe appears to be homogeneous. However, the 
possibility of anisotropic or inhomogeneous dark energy 



(or combined dark matter and energy sources) is begin- 
ning to be discussed in the literature [TtJl HT1 [T51 114) . 
Given our ignorance on the fundamental properties of 
these sources, there is no theoretically binding reason 
to assume, a priori, that no valuable new information 
would result from studying their interaction under inho- 
mogeneous conditions, at least in scales associated with 
structure formation and gravitational clustering. 

Alternative proposals to the dark energy paradigm also 
consider a full relativistic treatment of cosmological in- 
homogeneities [T5J [T3]. Numerous articles describe the 
possibility that cosmic acceleration (or at least part of 
it) could result from the presence of inhomogencites in 
photon trajectories in observations from high red-shift 
objects [TTl HE] . Following a more theoretical perspec- 
tive, various averaging formalisms |191 120] have exam- 
ined the occurrence of an effective acceleration from the 
so-called "back-reaction" , associated with non-local ef- 
fects that emerge as inhomogeneous sources are aver- 
aged (see [12 [53] for further discussion and [23J [25] for 
scalar averages in the spherical dust case). Also, an in- 
teresting connection between cosmic acceleration, back- 
reaction and non-trivial gradients of quasi-local energy 
has been proposed [26] . which could clarify important 
theoretical issues (the connection between the material 
presented in this paper and back-reaction issues is dis- 
cussed in [251 [2Z] ) ■ Possible observational consequences 
of these theoretical alternatives to dark energy are cer- 
tainly worth serious consideration [28] (see also [TTl HH] 
). Even if dark energy prevails over these alternative 
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proposals, there is no harm done (and perhaps valuable 
new information) in probing its behavior under inhomo- 
geneous conditions. 

Ideally, the study of inhomogeneous sources should be 
conducted on spacetimes not restricted by simplifying 
symmetries, but this would require fully 3-dimcnsional 
codes of high complexity. As a compromise between lin- 
ear perturbations and this type of "realistic" generality, 
we present in this article a class of inhomogeneous spher- 
ically symmetric models that can be well handled with 
relatively simple numerical methods. Although these 
models are idealized, they are non-trivial, exhibit non- 
linear behavior, and are general enough for testing a wide 
variety of physical assumptions on modern cosmological 
sources. 

By assuming a spherically symmetric Lemaitre- 
Tolman-Bondi metric, we derive in section |TT] a class 
of spacetimes ("LTB spacetimes") that generalize to 
nonzero pressure (and nonzero pressure gradients) the 
well know solutions for inhomogeneous dust associated 
with this metric [29j . The most general source for this 
metric (in the comoving frame) is a fluid with anisotropic 
stresses, which will end up being interpreted as fluctua- 
tions of the isotropic pressure. However, regardless of 
the interpretation, we remark that the existence of pres- 
sure anisotropy is far from a drawback or shortcoming, 
not only because it can be associated to a number of 
well motivated physical effects [3D], but because an in- 
homogenous source with only isotropic pressure is a far 
more idealized situation than one with anisotropic pres- 
sure. 

In section|TTT]we show that LTB spacetimes can be fully 
characterized by a set of local covariant scalars, hence 
their dynamics becomes completely determined by solv- 
ing the evolution equations for these scalars in the covari- 
ant fluid flow (or "1+3") representation |3T1 132"! 1331 134] . 
By looking at the relation between energy density and 
the Misner-Sharp quasi-local mass-energy invariant |35[ 
131)1 [3T1 [351 [35] , we introduce in section IV "quasi-local" 



(QL) scalar functions, which are LTB objects that satisfy 
FLRW dynamics. Expressing all local covariant scalars 
as perturbations of the QL scalars, we rewrite in section 
|V| the 1+3 fluid flow equations of section |III| as evolu- 
tion equations for the QL scalars and these perturba- 
tions, which on the basis of known criteria that define a 
perturbation formalism [3TJ [3SJ |33J [33] HJJ, are shown in 
section |VI| to be covariant, gauge-invariant, non-linear 
perturbations on a FLRW formal "background" given by 
the QL variables. These evolution equations are equiv- 
alent to ODE's in which radial dependency enters as a 
parameter (see Appendix C of [3D])). 

As shown in section |VII| once we choose an "equation 
of state" (EOS) between the QL energy density and pres- 
sure (which are "background" variables), the QL evolu- 
tion equations become determined. We point out that, 
by choosing such an EOS, correlation terms appear in 
the relations between local scalars, in a similar way as 
virial corrections to the ideal gas EOS appear in Newto- 



nian self-gravitating systems [H] H3J . These terms can 
be associated with the long range nature of gravity |44j . 
Since the fundamental physics of dark matter and dark 
energy sources is still unknown, we cannot rule out the 
possibility that this type of non-local effects could play 
an important dyamical role when these sources are stud- 
ied under non-trivial and non-linear inhomogeneity. 

In section |VIII| we apply the perturbation formalism 
to a fluid mixture, which can be interactive or with each 
component separately conserved (decoupled). This mix- 
ture description can be readily used to generalize to in- 
homogeneous conditions similar mixture models derived 
in a FRLW context in the many references quoted in 
[H] [7J. We show in this section that scalar fields 
can only be compatible with LTB spacetimes in mixture 
sources, playing the role of the homogeneous dark en- 
ergy component (coupled with inhomogeneous dark mat- 
ter [13]) 

Coordinate choices appropriate for the numerical treat- 
ment for the evolution equations of sections [V] and VIII 
are discussed in |IX[ while in section [X] we show that per- 
turbation functions are well defined as long as shell cross- 
ing singularities are absent. A possible (but not compul- 
sory) way in which LTB spacetimes can be embedded 
into a FLRW background is by smoothly matching an 
inhomogeneous LTB section with a section of a suitable 
FLRW spacetime. We provide in section XI the condi- 



tions for such a matching, leading to smooth and fully 
relativistic generalizations of the Newtonian "top hat" 
models that are widely used in the astrophysical litera- 
ture US]. 



We examine in section XII two examples of the applica- 



tion of the formalism presented in the previous sections. 
In the first example, we solve numerically the evolution 
equations for a smooth and fully relativistic Chaplygin 
gas "top hat" model, in which a black hole is formed in 
an inhomogeneous comoving section of the LTB Chaply- 
gin gas, smoothly embedded in an expanding Chaplygin 
gas FLRW universe. In the second example, we examine 
a mixture of interactive dark matter and dark energy of 
the type examined in [5J [6] UJ, but give the interaction 
a phenomcnological interpretation in terms of particle 
creation, with dark matter particles decaying into dark 
energy [Jj [46]. We show how the radial profiles of dark 
matter density evolve from initial clumps into deep voids, 
thus providing a toy model for a void formation scenario 
in the context of this type of interactions. In the ap- 
pendices we provide detailed information on how initial 
conditions and coordinate choices can be set up for the 
numerical integration of the evolution equations. 



II. GENERALIZED LTB SPACETIMES 

We will denote by "LTB spacetimes" all spheri- 
cally symmetric spacetimes that are solutions of Ein- 
stein's equations described by the Lemaitre-Tolman - 
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Bondi metric (in the comoving frame): 
R' 2 dr 2 



ds 1 



-c 2 dt 2 



1 



R 2 [dO l 



sin 2 9 1 



(1) 



where R = R(ct,r), R' = dR/dr and K = K(r). The 
metric is normally associated with a dust source, but 
the most general source compatible with LTB spacetimes 
is an anisotropic fluid [29 



T ab = fiu a u b +ph ab + U ab , 



(2) 



where u a — #Q,with n(ct,r), p(ct,r) being respectively 
the energy density and isotropic pressure, H ab is the 
anisotropic pressure tensor and h ab — u a u b + g ab is the 
induced metric of the hypersurfaces 3 T(t) orthogonal to 
u a , marked by constant t. 

The covariant objects associated with LTB spacetimes 
are the scalars /i, p, plus the expansion and Ricci scalar 
of the 3 T(t) (respectively 6 and 3 K) 



u a u b T ab , 



1 



P= 



h n ,T ab , 



6 = \7 a u a = 



2R R' 
~R + R' 



2{KR)' 
R 2 R' ' 



(3) 

(4) 
(5) 



together with the anisotropic pressure II a& , shear a ab and 
E ab Electric Weyl tensors 



IT 6 = 



i {a h b} 



V( a u a) 



&ab 



u c u d C ab 



\h ab h cd 


rpcd 


(6) 
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(8) 



where hi a h b J denotes symmetrization on a, b, R = 
u a V a R = dR/dct, V Q = h b a V b and C abcd is the Weyl 
curvature tensor. 

Each spacelike symmetric trace-less tensors like 
n ab , a ab , E ab in LTB spacetimes can be expressed in co- 
variant manner in terms of a single scalar function as 



n ab 


= VE ab , 


=► V = E ab Il ab , 


(9a) 


a ab 


= £S Q& , 


=> 1j = ^abCT , 


(9b) 


E ab 


= £E ab , 


=$> £ = E a bE ab , 


(9c) 



where the tensor E ab = h ab - 3x a X b and \ a = Vh rr 5° 
is the unit vector orthogonal to u a and to the orbits of 
SO(3) parametrized by (0,0). 

The scalars //, p, V relate to the metric functions by 
means of the field equations G ab — nT ab (with k = 



8ttG/c 4 ) for 0-(§: 

k /i R 2 R' = 
up R 2 R' = 
R' 



R(R 2 + K) 
1 r 



R 



R(R 2 + K) + 2R 2 R 
R 2 + K 2Y 



R 2 



Y 



(10a) 
(10b) 

(10c) 



while from Q and ([9]), we obtain for £ and E 

£ = 



R R[ 
R~R/ 



-V U 

2 6 ^ 



R 2 + K 
2R 2 ' 



(11a) 



(lib) 



It is evident that LTB spacetimes are completely charac- 
terized by the following set of covariant scalars (see [33] 
for comparison): 

{n, P , r, e, s, £ , 3 n}, (12) 

which we will denote as the "local" scalar representation. 



III. FIELD EQUATIONS VS. 

EQUATIONS 



FLUID FLOW 



In principle, LTB specetimes can be determined by 
solving the field equations after supplying an "equation of 
state" (EOS) linking //, p and V. This can be illustrated 
by looking at the best known LTB spacetime, a dust 
source with the simple EOS: p = V = 0, for which total 



energy density is rest-mass density: 
case, ( 10a I— ( 10c I reduce to 



k P c 2 R 2 R' = 2M', 



R 2 = 



2M 

~R 



-K, 



jU = per. In this 

(13) 
(14) 



where M — M (r) appears as an "integration constant" 
and the rest of the scalars O, E, 3 1Z, £ remain the same 
(with V = in ( lib I ) . The standard tactic in this case 



is to solve the Friedman equation (14) for a given choice 
of M(r) and K(r), and then to evaluate from this solu- 
tion all the derivatives of R, M, K in o rder to com pute 
the scalars p, 0, E, 3 K, £ from (Q, ((SJ, ( [Ha) , ( [TTb| and 
( 13 1. There is an extensive literature on how this is done 
(see [251 SB HH])i and the same tactic can work on sim- 
ple particular cases with non-zero pressure (see [551 US] ) 
and simple EOS relating p, p and V, but for a general 
LTB spacetime it is not possible (or very difficult) to find 
and solve an equation like (14) in order to evaluate radial 



gradients of R for computing the scalars in ( 12 1. 

An alternative approach to deal with LTB space- 
times follows from the evolution equations (or "fluid 
flow" equations), based on the covariant "1+3" de- 
composition of Ehlers, Ellis, Bruni, Dunsbury and Van 
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Elst [3TJ [351 [33J 131] ■ Instead of working on the relation satisfies: 
between ([l2l and the metric functions, in this approach 
the field and conservation equations V a T ab = for |l])- 
([2]) are transformed into an equivalent first order system 
of time evolution equations and constraints for ( |12[ ) . In 
particular, the spherically symmetric case is discussed in 

Em. 

Since the expansion and 3-dimensional Ricci scalars 
very often appear as 0/3 and 3 lZ/6, the notation be- 
comes greatly simplified if we define: 



H = 







K = 



6 ' 



(15) 



Bearing in mind ( 9a I— ( 9c I and (15), the "1+3" equations 



for LTB spacetimes become the following set of scalar 
evolution equations 



H 



6 



(li + 3 P ) -2E 2 



-3(/i + p)H-6E7>, 



E = -27fE + S 2 -£ 



[m + v - 



- -V 
2 ' 

2P]E 



(16a) 
(16b) 
(16c) 



3 (S + ^P) (H + E) 



together with the spacelike constraints 



(p-2P)'-6pi- = 0, 

3 S f = 0, 



(£ + «)'■ 



£' + 3£ 



R 



0, 



(16d) 

(17a) 
(17b) 
(17c) 



and the Friedman equation (or "Hamiltonian" con- 
straint) 



H 



- fi - JC + £ 2 



(18) 



The system ( 16a I— ( 18 1 still requires an equation of state 



(EOS) linking /1, p and V to become determined. How- 
ever, time and radial derivatives do not decouple, in gen- 
eral. Moreover, ( 12 1 is not the only scalar representation 



for LTB spacetimes. 



IV. THE QUASI-LOCAL (QL) SCALARS. 

The Misner-Sharp QL mass-energy function, A4, is a 
well known invariant defined for spherically symmetric 
spacetimes as [35 ^ [36 | l37 t [38 ] [39]: 



2M 
~R~ 



ee R 2 n e ^ M = 1 



{d r Rf (d R) 2 



9oo 



(19) 



where lZ ab d is the Riemann tensor. For LTB spacetimes 
with an anisotropic tensor like ([2]), M. defined by (19) 



2M' = kliR 2 R', 
2M = -n(p-2V)R 2 R 1 



(20a) 
(20b) 



whose integrability condition is equiv alent to th e co n- 
servation law Vf,T ab = (equations (16b I and (17a)). 
In fact, the scalars M and R are the main invariants 
in spherically symmetric sacetimes [351 EH], hence all 
objects constructed with R and M. and their covariant 
derivatives are necessarily covariant objects (see [36"1I3?] ). 
Comparing (10a I with (20a I suggests that M. can be 



equal to the integral of /i in (10a) as long as this integral 
can be performed and is finite. If we assume the exis- 
tence of a regular symmetry center at r = 0, then we can 



integrate both sides of ( |10a l along comoving domains of 
the hypersurfaces 3 T(i) to obtain 

2M_k _ k ^Zo^ R2R ' dx _ R 2 + K 
i? 3 _ 3 M * _ 3 r= r n R^R'dx ~ R 2 ' ( ' 

Jx—O 

which, for spherical symmetry, characterizes M. as the 
"effective" mass-energy in a comoving volume [37 ] l38 | l39] 
and the "ADM" mass in the appropriate asymptotic Xm\- 
its [351 [31 [371 [351 [3J [5U] . Notice that for the particular 
case of dust (p = V = 0) we have M. = M and the 
QL mass-energy is conserved along the 4-velocity flow: 
M = 0. When p, V are nonzero, then (211 generalizes 



the Friedman equation ( 14 1 . 

The function "/i*" that emerges from the inte- 



gral relation between \i and M. in (21) motivates us 



to define the following map acting on any scalar function: 

Definition. QL scalars. Let X be the set of all smooth 
integrable scalar functions in a comoving regular domain 
V = § 2 x £ C 3 T(i), where S 2 is the unit 2-sphere 
parametrized by (9, (j)) and £ = {x \ < x < r}, with 
x = marking a symmetry center. For any function 
A e X, a "dual" QL scalar function follows from the 
map J* : X — ► X such that 



^4* = J7*(^4) 



J^AR 2 R'dx 



(22) 



The QL scalar functions A* : T> — ► K depend on the 
upper integration limit r of £, and generalize to any 
scalar the QL mass-energy definition pi) constructed 
with ix. See [27] for a more detailed discussion and a 
covariant definition of the integral (22 1. 



QL scalars A* defined by ( 22 ) comply with the follow- 
ing useful properties: 



A A* 



1 



R 3 (r) J 



A'R 3 dx, 



(23a) 



(23b) 



i* = {A) t + 3{HA),-3H t A„ (23c) 
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where, in order to simplify notation, we have adopted the 
following conventions: 

A, = [JM)Y ^ JM) = (M*, 

AJ = [J.(A)\ ' + MA') = (A')*, 
(AB)* = J*(AB), 



. dx 



. dx, 



x=0 



(24) 



It is straightforward to verify that these properties are 
perfectly self-consistent. 

From @, ^ and the QL scalars dual to H, JC 

are 



h* = j*{H) = 

tC* = J*(/C) = 



R 

r: 

K 

r? 



/c* 



(25a) 
-2W*. (25b) 



Hence, by applying (22 1 to (10a I and (10b) we get 



k / pR z R'dx = = R 



(i? 2 + k) 



K I pR^R'dx = -p*R 4 = - 



(26a) 

U I R- + h' ) +2RR , 
(26b) 



so that, with the help from (25 1, these equations become 



identical to the Raychaudhuri equation of FLRW models, 
and its integral, the Friedman equation: 



H t . 



-Ul - - (/i* +3p*), 



Hi = — /i* — /C*, 



(27) 
(28) 



where H\ = [J*(H)] 2 ^ J*{H 2 ). These equations can be 
combined to yield the FLRW energy balance equation 



-3(/j* + p*)W* 



(29) 



Applying (221 to (10c), (Hal and ( lib I yields the re 



maining scalars, S, V and £ , as deviations of H, \x and p 
from their QL duals 



V ■■ 
£ = -'- 



■[H-7U], 



[p-p* 



(30a) 
(30b) 

(30c) 



while the QL mass energy function complies with 



2M = -(i*R 3 , 2M = -Kp*R 2 R. (31) 



V. EVOLUTION EQUATIONS FOR THE QL 
SCALARS. 

The local and QL scalars can be related by means of 
"relative deviations" 



5^ = 



^4* 



A = A* 



1 + 8^ 



(32) 



which allows us to eliminate /x, p, H in terms of their 
duals A* and the corresponding S^ A \ Hence, bearing in 
mind ( 30 1 , we have a complete scalar representation of 



LTB spacctimcs given by 



{H*, p*, P *, 5 {H \ 5<- p \ 8^}. 



(33) 



which is alternative to the local representation ( pL2[ ) . We 
will denote ( 33 1 the "QL scalar representation" . The 1+3 
system ( 16a I— ( 18 ) is the set of evolution and constraint 



equations for (12 1, we derive now the evolution and con- 



straint equations for the representation (33 1 



From ( 23b ) , the radial gradients of /x* , p* and H * can 



be given interms of the 8 functions 



H* 
H* 



3R' 
~R 



8^ H \ 



A** 



3K 
~R~ 



8^\ 



p* 



3R' 
~R 



(34) 

while ( 27 1 and ( 29 1 are evolution equations for ji* and 
H*. Hence, the evolution equations for 5^*) and 8^ 
follow from the consistency condition 



[A*'}-= A 



(35) 



applied to (27 1, (29 1 and (34 1 for A = H*, (i*. The result 



is the following set of autonomous evolution equations for 
the QL representation ( 33 1 : 



= — 3 [ 1 + w] /i* H*, 

= -Hi- j. [1 + 3™] p., 



(36a) 
(36b) 



8^ = m 



(>)_£(?)) w - (l + w + tffAsM 

(36c) 

7U (l + SW) 8M 



6H* 



where 



p* 
ft* 



(36d) 



(37) 



The spacelike constraints associated with these evolution 
equations are simply the spatial gradients p4| , while the 
Friedman equation (or Hamiltonian constraint) is (28 1. 
Equations (36a)-(36dl become fully determined once an 
"equation of state" that fixes p*, 8^' as functions of 
/i*, is selected (see section VII). 
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With the help of (30 1 and (32 1, it is straightfor- 



wards to prove that the evolution equations ( 36 I and the 



constraints (28 1 and (34 1 are wholly equivalent to the 
1+3 evolution equations (16 1 and their constraints (17 1. 



Hence, given an "equation of state" , they completely de- 
termine the dynamics of LTB spacetimes. 



VI. 



A NON-LINEAR PERTURBATION 
SCHEME 



Th e definition (32 1 and the evolution equations (36al- 
( 36d I suggest that the role of the functions 5^ can be 



rigorously characterized as spherical perturbations on a 
formal FLRW "background" state given by the A* . The 
perturbation formalisms developed by Ellis, Bruni, Duns- 
bury and van Elst [3U 02 [33j E2] and Bardeen [41] define 
a perturbation scheme in terms of a suitable map $ be- 
tween objects A in S (a FLRW model) and objects A in S 
(the "perturbed" lumpy model). The perturbation then 
compares A with which are objects in S that define 

a formal FLRW "background" . The application of this 
formalism to the case when the lumpy model is an LTB 
spacetime simplifies considerably because (i) both classes 
of spacetimes (LTB and FLRW) can be completely de- 
scribed by covariant scalars, and (ii) both are given in the 
same normal comoving coordinate (or frame) representa- 
tion. Hence, 4> can be a map between covariant scalars 
in S and covariant scalars S, and practically all problems 
due to gauge and coordinate freedom disappear. 

The map $ can be defined as follows: let X and X be, 
respectively, the sets of smooth integrable scalar func- 
tions in S and S, then for all covariant FLRW scalars 
A E X (we denote FLRW objects with an over-bar), $ 
is the map 



<$> :X->X, $(A) = J*(A) = A* e X, 



(38) 



which characterizes QL scalars {p*, p*, H*, 3 7Z*} (LTB 
objects satisfying FLRW dynamics) as the "background" 
model in LTB spacetimes, while the perturbations 5^- A ' 
provide their comparison with covariant scalars A 



§ (A) 



A-$(A) 



(39) 



Following Dunsbury, Ellis and Bruni [3TJ 152] , a per- 
turbation scheme on FLRW cosmologies is covariant if S 
is described by the 1+3 fluid flow variables, as in (16) 
[3H |33|331 131]. Although our description of LTB space- 
times is not based on these scalars (the representation 
( 12 1 ) , it is still covariant because {p*, p* and ©*}, are 



themselves covariant scalars by virtue of their connec- 
tion with the invariants Ai, R and their derivatives in 



( [25| and (311 (see [37])- Hence, the formalism associated 
with (38) and (39 1 is covariant. 



Also, by virtue of the Stewart-Walker gauge invariance 
lemma |31j . all covariant objects in S that would vanish 
in the background S (a FLRW cosmology in this case) are 



gauge invariant (GI), to all orders, and also in the usual 
sense (as in [41]). The background variables /i*, p*, ©* 
do not vanish for S, hence they are "zero order" GI vari- 
ables to all orders. The quantities in LTB spacetimes that 
vanish for a FLRW cosmology are the tensors n ab , a ah 
and E ab , given by (9| in terms of the scalar functions 
V, £ and £ in (flOc l-([Tl| . But from Q, these func- 
tions are basically the fluctuations p — /u*, p — p* and 
— 6*. Hence, from (32 1, (34 1 and (39 1, the perturba- 
tion variables 6^ , S ^ and S ( 1 , as well as the gradients 
fi'^ , p'* and 0* , are all "first order" quantities that are 
GI to all orders. Therefore, LTB spacetimes in the QL 
scalar representation {A*, 5^} can be rigorously char- 
acterized as spherical, non-linear GIC perturbations on 
a FLRW background. In the linear limit \S^ \ <C 1 these 
perturbations reduce to spherical perturbations in the 
long wavelength approximation and in the synchronous 
gauge [27] . 



VII. EQUATIONS OF STATE. 



The system (36 1 is still undetermined, as there are no 
evolution equations for p* and S^ p \ As in any perturba- 
tive approach, we need to impose an "equation of state" 
(EOS) between the background variables, which are in 
this case p* and /x* . Such a n EO S determines the "back- 
ground subsystem" (36a)-(36bl, as well as the relation 



between the perturbations 5 {M> and 5^ p \ 

Consider a commonly used barotropic equation of the 
form 



P* = P*(p*)- 



(40) 



From (34 1, we obtain the corresponding relation between 



fluctuations and perturbations of p and p 



P-P* 



dp* 
dp* 
dlnp* 



[p, — /i* 



d\n 



(41a) 
(41b) 



Notice that, in general, the EOS (40 1 will not be satisfied 



by p and /1. But this is expected, since in a perturbation 
scheme fluctuations or perturbations do not satisfy, in 
general, the EOS of the background (not even in linear 
perturbations). The exception is the simple particular 
case of ((40 1 



p* 



w (i* 



(42) 



where wq — ino(i), and in particular wq can be a constant 
(so that w = is dust). The EOS (E2} leads to 



p - p* = w [ p - [1* 



(43) 



and so p = wq At, even if wq — wo(t) because the map J* 
in ( 22 1 involves integrals of scalars along hypersurfaces 
of constant t, hence J*(wo(t), p) = wo(t) J*(p) for all t. 
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It could be argued that, because the perturbations are 
"exact" and non-linear (i.e. S^ A ' can be "large"), even if 
a given EOS ( 40 1 is physically reasonable we could have 



a wholly unphysical relation between p and \i through 
(41a I. While this possibility cannot be ruled out, it is 



something that must be tested by looking at the solutions 
of (36 1 for different EOS. Evidently, the plausibility of the 



formalism must be verified and judged for different EOS 
according to the predictions of the solutions of ( 36 1 . 



However, it would be wholly incorrect to disregard the 
formalism from the outset on the grounds that local vari- 
ables /i and p do not satisfy a given "familiar" EOS (say, 
the Chaplygin gas) used in a FLRW framework. Not 
only we are treating /x and p as non-linear perturbations 
(not expected to satisfy the EOS of the background), 
but also relations that could be valid in a FLRW con- 
text (where local and QL variables are identical) need 
not hold automatically in an inhomogenous non-linear 
context. Demanding that the local EOS must be con- 
sidered as the only physically meaningful EOS in an in- 
homogeneous context is only correct for hydrodynamical 
sources characterized by short range interactions or for 
self-gravitating systems in a kinetic regime [44] . There 
is no reason to impose such condition in self-gravitating 
systems dominated by long range interactions. 

As shown in Newtonian self-gravitating systems in spe- 
cific scales [12J [53], the correlation effects of the long 
range nature of gravity appear as a first order virial cor- 
rection to the classical ideal gas EOS linking /i and p. By 
looking at (41a I, it is evident that an EOS like (40) in- 



troduces similar correction terms to the relation between 
local variables /j, and p (effects of this type in the classical 
ideal gas are discussed in [27] ). These corrections could 
convey important non-local information in the study of 
self-gravitating systems, as for example the so-called 
"back-reaction" , which can be related to dynamical and 
observational effects of QL energy [22] [23] [26]. Non- 
local correlation terms could also be important in look- 
ing at self-gravitating systems with relativistic dark en- 
ergy sources, as there is no reason to assume that these 
sources (whose fundamental nature is still unknown) are 
of a hydrodynamical nature even under inhomogencous 
and non-linear conditions. 



VIII. MIXTURES OF DARK MATTER AND 
DARK ENERGY. 

A. General mixtures 



The resulting models provide an inhomogeneous gener- 
alization to known DM-DE mixtures in a FLRW con- 
text [5] El [U [S]. To examine a DM-DE mixture we as- 
sume a decomposition of the "total" energy-momentum 
tensor as a sum T ab = T ab + T ab , which implies that 

DM DE ' 1 

energy density and pressures split as 

M = Mdm +M de > ( 44a ) 
P = P DM +P de , (44b) 
n ab = VE ab = [P nM + V^]E a \ (44c) 

where S ab was defined in In order to simplify nota- 
tion we introduce: 



P ~ A*DM ! 

<7 = fe, 
= p* + q* 



V = Pi>u, 

K = Pde > 

p* = (fi* 



while, from (30b I, V splits as in (44) 



V 



(45) 



(46) 



The total conservation law V T ab = can be split in 
conservation laws for the individual tensors 



v T ab = -a = _y rpa 

dm J ° ni 



6 

DE 1 



(47) 



where j a is the interaction current (or coupling) that 
characterizes an interactive mixture, so that if j a = 
the mixture is non-interactive (decoupled) . We take this 
current as a vector parallel to the 4-velocity, so that 
j a = Ju a and h ca j a = hold. The spatially projected 
conservation equation h ac V T ab = remains as it is, but 
the projection along u a becomes 



u«V 6 T°* = J = -u a V b TZ 



b 

DE ' 



(48) 



which, from (16b), (|44|) and (451, takes the form 



p + 3(p + (p)H-3(H-K,)(<p-<p.) = J, (49a) 
q + 3(q + ir)H - 3(H - H,)(tt - 7r*) = - J, (49b) 



where we used (30) and B6|) to eliminate V and S in 



terms of the fluctuations of <p, it and TL. With the help 
of (23c I and the relation between local and quasi-local 



variables (32 1, and considering that 



(50) 



The formalism that we have presented can be extended 
to describe sources like fluid mixtures, for example, a 
mixture of dark matter (DM) and dark energy (DE). 



we can transform ( 36 1 into the following set of evolution 



equations for an interactive mixture in the QL scalar rep- 
resentation: 
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n,, 
p* 

q* 
S(p) 



-Hi - - [p*(l + 3w p ) +q t (l + 3w q )} 

-3/9* [l + w p ] Hi, + J*, 
-3g» [1 + w 9 ] H* — J*, 



■J 



\w p (<5 (p) - 5 M ) -(l + w p + frhsM] - — (6^ - S (J) ), 

P* 

= m4w q (s {q) -s {7T) )-(i + w q + s^)s^}-—(s^-s^ J '>), 

SM = -HjW(l + 6W) + -^- 



(51a) 

(51b) 
(51c) 

(51d) 
(51e) 

* 

(51f) 



where 



1 

P 



w u 



7T 
<l' 



(52) 



The evolution equations for a non-interactive mixture 
follow by setting J* = <5( J ) = 0. An important partic- 
ular case follows by considering DM as cold dark mat- 
ter (CDM), hence = T»£ M takes the form of dust 
(ip = ip* =0). We examine numerically an interactive 
mixture in section IXIIBI 



B. Scalar fields. 

The scalar field ( "quintessense" ) is a very popular 
model for dark energy sources. Unfortunately, the 
energy-momentum tensor of an inhomogeneous scalar 
field is incompatible with the LTB metric. Hence, the 
only way in which we can model a scalar field with LTB 
spacetimes is by considering it as homogeneous DE, while 
DM is then an inhomogeneous fluid that can be (if de- 
sired) CDM dust, but this is not a necessary condition, 
though if we take DM as CDM, then T ab takes the form 
of a perfect fluid {V — 0). 

The particular ca se of a mixture with homogeneous 
DE follows from {FT} with 5^ = <jW = 0, so that q = q« 
and 7r = 7T* , and it implies the following constraint on 
the interaction term </* 



Ji = 3q,(l+w q )K 



=S> J,S {J) = 3q(l+w q )Hj( n) , 

(53) 

so that uncoupled mixtures are not possible (J* = im- 
plies a FLRW spacetime 5^ = 0). In particular, we can 
always assume a scalar field form with arbitrary potential 
for the DE homogeneous fluid: 



v(<f>), 



(54) 



with <fi = </>(t), which transforms (51c I into the Klein - 
Gordon equation 



37i* 



1 



5W \ 
§W J 



dV 



= 0, 



(55) 



The mixture examined in |13j is the particular spatially 
flat case of such a mixture, in which the interaction term 
J* was obtained from the Lagrangian of a scalar-tensor 
theory. 



IX. A FLRW— LIKE LTB METRIC. 

It is useful for the numerical solution of the evolution 
equations to re-parametrize the LTB metric ([TJ so that 
it looks as close as possible to a FLRW line element. 
This metric is invariant under arbitrarily re-scaling of 
the radial coordinate r = r(f), hence we will use this 
gauge freedom to chose this coordinate so that 



R i = R(ct i ,r) = h- 1 r, 



(56) 



where h s 1 is an arbitrary characteristic length scale. 
Considering ( 56 1 , we introduce the following dimension- 



less form for the metric functions R and R': 

e = ^ = ^ l = i, 



r 



R R 

Ri^K 1 
R'/R 
Rj/R~ 



1 



rt' 



r, 



(57a) 
(57b) 



where the subindex j will henceforth denote evaluation 
at a fiducial hypersurface 3 7^ marked by t = i,-, which 
can be considered as an initial surface for treating the 
integration of ( 36 1 and pTj ) as a well posed initial value 
problem. Notice that [vC]7 = [Ai]„ holds in general. We 
can now write ( 34 1 for any A* as 



A.. r 



(58) 



The parametrization given by ( 56 1 and ( 57 ) transforms 
the LTB metric into the following form that resembles 
a FLRW line element 



ds 2 = -c 2 dt 2 + L 2 



T 2 dr 2 
1 - N;r 2 



{d9 2 



(59) 
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where 



6 hi 



hi 



(60) 



In terms of £, the QL scalars in (15) and (25) take the 
following FLRW-like forms 



e* 
T 

6 



I _ L 
I ~ V 
[fC*]i _ [k* 



L 2 



(61) 
(62) 



which highlights the role of £ (or L) as a position de- 
pendent scale factor, as opposed to R, whose geometric 
meaning is the "area distance" or curvature radius of the 
orbits of SO(3). 

The radial coordinate choice associated with R4 in ( [56] ) 
is appropriate when we have a symmetry center marked 
by r = 0, but then it is perfectly valid to set the radial 
coordinate so that 



(63) 



where f(r) is any smooth positive and monotonously in- 
creasing function so that /(0) = 0. Such choice was made 
in the mixture example of section XII B (see Appendix 
D). Other choices of f(r) depend on the topological class 
of the hypersurfaces 3 T, for example, if they have S 3 
topology (like a "closed" FLTW cosmology) , then a good 
choice for f(r) would be a sine-like function with a sec- 
ond zero at the second symmetry center. For wormholc 



topologies, S 



or § x S , the function / should not 



have zeroes but must have at least a 
where /' = 0. See 001 El HE] ■ 



'turning value" 



with £(t, r) playing the role of a position-dependent 
dimension-less FLRW scale factor a(t)/a(t ) and [£*]$ 
is a radially dependent "curvature index" . Once an EOS 
like (40 1 is given, equation (64b I can be formally inte- 



grated leading to 



= fx*(£, [fJ,*]i), 



(65) 



where = fj.*(ti,r) is the initial QL energy density 

(the "arbitrary constant"). Once (651 is known, the 
Friedman equation (64c I can be, in principle (and for 
simple enough EOS), integrated or reduced to quadra- 
tures 09"] . 

While considering the integration of the FLRW sub- 



system (64 1 (or using known solutions of it) is useful in 



itself, given the vast amount of known FLRW solutions 
[H El HI El El E] , the resulting integrals are far from a solu- 
tion to ((36l and d5ll. As opposed to the FLRW case, ((64 1 



are PDE's (not ODE's) and so initial values arc functions 
of r (not constants), hence there is a non-trivial radial 
dependence of all variables that must be still determined 
by solving the evolution equations for the &( A \ which re- 
late to radial gradients of the A* through p4| and convey 
the departure from homogeneity. Still, the solutions to 
this subsystem are useful to discuss and understand a ba- 
sic regularity condition: absence of shell crossings, which 
is also important in LTB dust solutions [29l l40l 1371 148] . 

Since /x* in (651 is, in general, independent of T, we 



can use use ( 57b I and to eliminate £' in terms of T and 



applying (58 1 to this form of n*(£, [/x*]j) we obtain the 



general relation 



1 

3T 



(r-i) 



91n/i* 



35, 



(m) 



d In /x* 
dln[/x«] 



(66) 



X. 



THE "HOMOGENEOUS SUBSYSTEM" AND 
REGULARITY CONDITIONS. 



The evolution equations ( |36a l and (36b I for the back- 
ground QL variables /x*, 7i» are formally identical to 
FLRW evolution equations (like the evolution equations 
for q x , in (51 )). This suggests exploring the possi- 



bility of integrating separately these "homogeneous sub- 
systems", or even using known FLRW solutions, to de- 
termine analytically the background variables in ( 36 1 and 

dSIi. ' ' 



The metric form (|59|), with (|61|) and (62 1, allows us 

-(36b I and the 



to re-write the evolution equations (36a 



Friedman equation ( 28 1 in the following familiar FLRW 
form: 



[/x* + 3p» 



i 

i>2 



A«* +P* 



(64a) 

(64b) 
(64c) 



It is evident from ( 66 1 that 5^' and fi diverge for finite /x* 
if r — > for I > 0. Therefore, a necessary and sufficient 
condition to prevent this unphysical feature is to demand 
that initial conditions are selected, so that the following 
constraint: 



r(t,r) >o, 



(67) 



holds for all (t, r) for which £(t, r) > (the range of 
physical evolution) . This is equivalent to demanding ab- 
sence of shell-crossings associated with r = occurring 
in the physical evolution range 001 HZ1 HE] ■ This regular- 
ity condition applies to all scalars, such as and 6^ n \ 
and in general prevents an unphysical evolution in which 
perturbations 5^ A ' diverge, hence making local scalars A 
diverge, all of it while their associated QL variables A* 
remain finite. 

Initial conditions that fulfill ( [67| can be given ana- 
lytically for dust sources 03 148} 151] . In general, this 
regularity condition must be verified numerically for any 
choice of initial conditions. 
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XI. SMOOTH AND RELATIVISTIC "TOP HAT" 
MODELS. 

When assuming the existence of a symmetry center, a 
type of model that is often considered is that of a single 
spherical clump or void region that could be, somehow, 
embedded in an asymptotic homogeneous FLRW back- 
ground. The non-linear perturbation scheme that we 
have presented is ideally suited for this type of scenar- 
ios, and it does not restrict the different ways in which 
boundary conditions can be set for this asymptotic blend- 
ing into a cosmic background. 

One possibility for setting boundary conditions is sim- 
ply to perform a smooth match, at a fixed comoving ra- 
dius r = rb , between a region of an LTB spacetime (con- 
taining a center) and a suitable or "equivalent" FLRW 
spacetime (the precise meaning of "equivalent" will be- 
come clear further ahead). The resulting compound 
spacetime is formed by a comoving LTB region (a clump 
or void) in the range < r < rt, smoothy matched to a 
homogeneus FLRW section that extends for r > r^. This 
construction generalizes to smooth and relativistic con- 
ditions the so called Newtonian "top hat" or "spherical 
collapse models" that are widely used in astrophysical 
literature [45], and that tend to overlook the details of 
the embedding of the overdensity or void into the FLRW 
background. 

The section of a FLRW spacetime to be matched with 
a LTB comoving region is described by the metric 



ds 2 = 



2 dt 2 



dr 2 



(de 2 



1 — kor 2 

(68) 

where ko — 0,±1, and r > r b - This FLRW section is 
characterized by the covariant scalars {p, p, 7i, 1C}, with 
TL = a/a and K. — ko/a 2 (FLRW scalars are denoted by 
an overbar). The metric of the LTB region is just (59) 
with < r < ?v Since both metrics are already given 
in the same normal comoving coordinates, the conditions 
for a smooth matching between them simplify consider- 
ably. 

The standard matching conditions between two space- 
times are those of Darm ois 1471151] . For the LTB-FLRW 
case, with the metrics (68 1 and (59 1, these conditions 



imply |47j continuity at r b of their metric function 
and its proper time derivative (which is coordinate time) . 
From (57 1, (59 1, (61 1, (62) and ( [68] ), these requirements 
lead to 



pressure in the radial direction, T r 
sides (69), we have from ( 30b[ ) 



p — TP , hence be- 



\p-2V} h = \p*} h = p(t). 



(70) 



Since in the FLRW region QL and local scalars are 
identical: p = p*, in the LTB region we must have 
[p — p*]b — > as r — > rb for all t. Hence, (|30b|, p4|) and 



( 70 1 clearly imply 



V h = 0, p'*(t,r b ) = 0, \p-p*]b= [<5 (p) ]b = 0, 

(71) 

so that Q takes the perfect fluid form at r = r b - 

Since Darmois matching conditions require continu- 
ity of the QL scalars {/x*, j?*, /C*} and the FLRW 
scalars {p, p, K.}, the "equivalent" FLRW spacetime to 
be smoothly matched to a given LTB region can be pre- 
cisely defined as that in which p and p satisfy the same 
EOS that p* and zx* satisfy in the LTB region. Also, from 
(62 1, ( 69a I and ( 69d I , QL spatial curvature must be con- 
timious: |/C*] b = k /a 2 (t). Hence, the spatial curvature 
of the whole matched FLRW region is determined by the 
sign of [/C*] b . Notice that /C* changing sign for values 
r < rb in the LTB region has no effect on the matching. 

While Darmois matching conditions do not require 
continuity at r = of the local energy density ix, nor of 
the metric component g rr , it is highly desirable to avoid 
discontinuity of these quantities. Hence, we supplement 
these conditions by further demanding that 



Mb 



o, 
l. 



^(t,r h )=Q, (72a) 
(72b) 
(72c) 



The complete set of matching and regularity conditions 



( 69 1— ( 72 ) lead to a well defined and self -consistent model 



of a spherical comoving LTB region that is smoothly 
matched to a appropriate FLRW background. These con- 
ditions imply: 



§(p) 



5 (H) 



= 0, 



(73) 



and consequently their time derivatives vanish as well. 
The system (36) reduces at r = r b (and for r > rb) 



to the FLRW equations ( 36a I— ( 36b I , while the mixture 



equations (51 1 reduce in the same ranges to the FLRW 
equations ( 51a H 51c). 



i k *]i Ob) 



a(t), 
H(t), 
Pit), 
k , 



(69a) 
(69b) 
(69c) 
(69d) 



where the subindex b indicates evaluation at r = tv 

For a tensor T ab with anisotropic stresses like Q , Dar- 
mois conditions also imply |52] the continuity of the total 



XII. APPLICATIONS 

We examine in this section two numerical applications 
of the perturbation formalism: black hole formation in 
a Chaplygin gas universe and evolution of voids in a 
mixture of cold dark matter (CDM) and dark energy 
(DE). These examples are idealized because their pur- 
pose is to illustrate how the formalism can work in prac- 
tice. More "realistic" and comprehensive studies, and/or 
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models with a careful fitting to obervations, would re- 
quire separate articles in their own merit. 

We still need for the numerical solution of the evolu- 
tion equations ( 36 I and ( |5T| to re-cast the QL variables 
in dimension-less form (the perturbations 5^ <5 (p) , 5 {H) 
are already dimensionlcss) . Dimensionless quantities can 
be constructed by means of the freed om afforded by the 
arbitrary length scale h^ 1 defined in (56 1: 



/'* 



p* 



7U 

h. 



(74) 



which can also be done with the QL variables in (51 1 



The evolution time parameter, ct, must then be replaced 
by the dimensionlcss time 



h s c(t -ti). 



(75) 



so that t = corresponds to 3 %. The initial conditions 
and other technical issues connected to the numeric in- 
tegration of evolution equations in this section are pre- 
sented and discussed in the Appendices. 



A. Gravitational collapse in a Chaplygin gas 
universe 

The Chaplygin gas can be studied as an LTB spacetime 
in a fully relativistic inhomogeneous context. So far, this 
source has been examined under inhomogeneous condi- 
tions only in a Newtonian context |10j . as non-linear but 
relativistic approximations [llj . or as a static object [12] . 

The LTB spacetime associated with the "standard" 
Chaplygin gas EOS [31 H] is that characterized by the 
following particular case of (40 I 



a 



9 ■ 



(76) 



where a is a constant. Following (41), the relation be 



tween pressure and density fluctuations and perturba- 
tions is then 



p-p* = -zbi -(!»), 



(77a) 
(77b) 



sp ecialized for 

the EOS ([76} and with $0») given by~f77b). Since p = 
p*[l + 5( p> \ and fj, = + 5^], equations (77) imply 
the relation 



P 



■ a [l 
A* 



(78) 



which can also be interpreted as p being expressible in 
terms of the familiar (local) EOS, with a first order virial 
correction due to squared fluctuations that convey the 
effect of long range interactions. This type of relation 
is qualitatively analogous to that between the state vari- 
ables of the ideal gas under a self-gravitating regime in 



Newtonian systems [HJ U3J (see also [57] ) . As long as we 
ignore the fundamental physics of the Chaplygin gas, us- 
ing the EOS ( 76 1 and having the local variables p, /i given 



by (78 1 in an inhomogeneous self-gravitating context, is 



a reasonable working hypothesis. 

The solutions to the homogeneous subsystem (64 1 
leads to the following forms of and p* as functions 
of £ and [/z*]. 



A** = 



P* 



t 3 



1/2 



{a(*»-l) + M} 1/a ' 



(79a) 
(79b) 



These equations clearly illustrate the well known behav- 
ior of the Chaplygin gas as dust-like (p* ss 0) for earlier 
times and as cosmological constant for latter times: 



{[/yrj - a} 

i 3 



1/2 



d 3 



fa, i<l 



1/2 ' 



(80) 
(81) 



These asymptotic limits suggest considering the numer- 
ator of /i* in (80) as a rest-mass density, while (81 1 sug- 



gest considering ^fa as a cosmologic al c onstant. Using 
the length scale hj 1 introduced in (56 1, we define the 
following dimensionless parameters 



A = 



3/i2 



1/2 



(82) 



which lead to the dimension-less form for the initial value 



functions (see equations (99l-(100l in Appendix C). 

Ap plying the Chaplygin gas EOS to (66) and using 
( 79a I , we get the following expression for d^' 



5M = 



L-l) a + (l + ^ ) )[/z 

W 6 -i)«+W}r 



i. 



(83) 



which shows that F 
p diverge for finite jx 



> for I > implies that [i and 
and (which do not depend on 
F, as can be seen from (79 1). Hence, the condition (67 1 



to avoid shell crossings must be imposed. In principle, it 
could be possible to use (79a I for solving the Friedman 



equation ( 64c ) , and use the solution to find T (from 57b )) , 
and then S (fl> from (83), however it is far easier (and 
more intuitive) to solve the numeric system ( |36[ ) for the 
Chaplygin gas EOS. 

We will consider a "top hat" model, as discussed in sec- 
tion XI formed by an inhomogeneous Chaplygin gas co- 
moving region, matched at a comoving boundary r = 
with a Chaplygin gas FLRW background (its equivalent 
FLRW model). Assuming the initial inhomogeneity as 
that of an early universe small over-density with posi- 
tive but small scalar curvature, expecting that comoving 
layers in the overdensity will expand and then re-collapse 
while the FLRW background keeps expanding. 
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Since initial conditions depend on the parameters 
[m*]j, A and and we are assuming that t = U 

(t = 0) corresponds to nearly homogeneous early uni- 
verse initial conditions, then we can take h s — Hi, where 
Hi is a FLRW Hubble length at t = ti. Hence, by iden- 
tifying A = \J~A and bearing in mind that in the present 
era t = t$ we have 



kA 
3H$ 




0.7, 



(84) 



thus, A ~ 10~ 3 <C 1 since Hi ^> Hq. On the other hand, 
nearly homogeneous and spatialy flat initial conditions 
imply [ht]i ~ and [m*^ ss 1. The Chaplygin LTB plus 
FLRW system must comply with the matching conditions 
d69j-(|72]) to a spatially flat FLRW background (k = 0). 



Considering = 1, the overdensity extends in the range 
< r < 1 and the Chaplygin gas spatially flat FLRW 
background extends along r > 1. 




FIG. 2: Density contrast. The figure depicts the ratio 
fi,*(T,r)/p,(T) for a Chaplygin gas overdensity (0 < r < r h = 
1) smoothly matched to a Chaplygin gas FLRW universe 
(r > 1) with density p,. The matching interface r — = 1 is 
displayed as a white strip. The density contrast diverges as 
inner layers of the overdensity collapse into a black hole, while 
external layers blend into the cosmic background as r — > 1. 




FIG. 1: Hubble expansion and collapse into a black 
hole. The figure displays the function arctanTi* (r, r) for r > 
corresponding to a Chaplygin gas overdensity (0 < r < 1) 
smoothly matched to a Chaplygin gas universe at r = rb = 1 
(white strip). The dimension-less time r is defined in ( |75[ ) 
so that t = is the initial hypersurface. The function Ti, 
diverges at the initial big bang singularity. For inner layers it 
becomes negative and goes into Ti, — > — oo as these layers col- 
lapse into a black hole. For external layers blending into the 
cosmic background TL, —> yX, indicating an asymptotically 
de Sitter behavior. 



As expected on intuitive grounds, inner layers in the 
inhomogeneous overdensity bounce and collapse, as the 
FLRW region keeps expanding. This is shown in the 
graphs of figures [T] and [2] The kinematic evolution of the 
overdensity follows from the behavior of the QL Hubble 
scalar, — 0*/3 = £/£. As shown in figure[l] this func- 
tion is initially infinite at early times for all r (overden- 
sity and background), indicating an infinite expansion 
rate from an initial big bang singularity at £ = 0. As 
r increases and £ grows, the function H* monotonously 
deceases in the background and external layers of the 
overdensity, but in inner layers we have 7i* — > 0, which 




FIG. 3: Energy density fluctuations. The figure displays 
the relative fluctuation, or perturbations, 5^ (r, r) for r > 
corresponding to a Chaplygin overdensity smoothly matched 
to a Chaplygin gas universe at r = rb = 1 (white strip). This 
function is negative (overdensity) and close to zero near the 
center and the matching interface (small radial gradient). It is 
significantly different from zero in the areas where the radial 
gradient is large and <5' M - ) — + — 1 as inner layers collapse into 
a black hole and at the initial singularity for all layers (not 
shown) . 



means a zero of £ indicating that inner comoving layers 
near r — reach maximal expansion. See equations (61 )- 
( 64 1 and ( 79 1 . For later times (around r ~ 40 — 50j~wc 
— oo for these layers, indicating an infinite 



have 7i* - 

contracting rate associated with a re-collapsing singular- 
ity: curvature radius and scale factor, R and £, tend to 



zero, while /i* diverges (from (79al) 
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FIG. 4: Fluctuations of the Hubble scalar H. Fluctua- 
tions of Tt given by TL - TL, = H*8 (n) . As in figure |] it is 
close to zero near the center and the matching interface and 
it is significantly different from zero in the areas where the 
radial gradient is large. Notice that AH. = — E, where E is 



the shear scalar function in ( 9b I and ( 30a l 



The diverging of /i* for these layers is shown in fig- 
ure [2j which displays the ratio of the density /i* in the 
inhomogeneous region to the cosmic background density 
£l (the so-called "density contrast" in Newtonian astro- 
physical applications jH]). Notice that diverges 
at the same locus in the inner part of the overdensity 
where Ti* — > — oo, while external layers of the overden- 
sity smoothly blend to the background value: p,/p — > 1. 

The evolution depicted by figures [T] and [2] is a clear 
signal of a "big crunch" singularity that forms a local 
black hole where Ti, — > — oo and fi — ► oo. On the other 
hand, for the external layers of the overdensity TL* re- 
mains always positive and we have 7i 2 — > A ~ 0.001 and 
/i* — > A 2 ~ 10 -6 as t — > oo as they blend into the cos- 
mic background at r = 1. This asymptotic behavior is 
consistent with the expanding Chaplygin gas taking in 
this future asymptotic limit the form of a cosmological 
constant. 

Relative fluctuations of the energy density, are 
shown for r > in figure [3] This function is identically 
zero for r > 1 since radial gradients are zero at the cosmic 
background. It is negative but very close to zero in areas 
with negligible radial gradients (near the center r = 
and near r = 1) . It is significantly different from zero for 
intermediate values of r where radial gradients are not 
negligible, while 5^ — > —1 as the internal layers collapse 
into a local black hole (compare the locus of 5^ — > — 1 
with that of TL, — > — oo and fi* — > oo). Since fi = + 
5^'], this behavior apparently implies that p — > as /i* 
diverges, which would be wholly unphysical. How ever , 
— * —1 does not imply p — ► 0. As shown by (83 1, 
— > — 1 if T — > oo. In fact, numerical numeric tests 
show that r — > oo as I — > (big bang and re-collapse 
singularities), but — > in these limits. As I — *• 
both and /x* diverge because the product + iW] 



diverges even if 1 + 5^ — > 0. 

It is important to mention that pressure anis otrop y can 
be directly related to 5^*) by means of (30b I and (77 1, 
leading to the ratio 



V 



2(1-JM) ! 



(85) 



hence, pressure is nearly isotropic when radial gradients 
are negligible and reaches its maximal anisotropic ratio 
V/p = l/4 as -> -1 in the singularities, which 

can then be characterized as anisotropic (as well as non- 
simultaneous) . 

The relative fluctuations of the Hubble scalar — 
(TL — TL,)/TL, will necessarily diverge for inner layers at 
their maximal expansion (TL, = 0). Hence, we plotted in- 
stead in figurejijthe fluctuations TL—TL, = H*5^ = — S. 
As shown by this figure, TL — TL, = for r > 1 (fluctu- 
ations vanish at the cosmic background) and is negative 
but small for most of the range of r and r ("clumps" 
of TL), but as TL, (see figure [TJ TL — TL, — > — oo at the 
collapsing singularity. 



B. Void formation in an interactive CDM-DE 
mixture 

A phenomenological interactive CDM-DE mixture in 
a FLRW background was examined in the three articles 
of reference [B] and reviewed more recently in [7j. This 
model assumes that the CDM component is dust and the 
DE component is a fluid satisfying a simple linear EOS 
of the type (42 1. It corresponds to the particular case of 



( 45 1 given by 



<p, — S v = 0, 
tt, = wq q*, 



w p = 0, 



$("•)=$(«), w q = w , (86) 



where wo is a constant in the range — 1 < wq < —1/3. 
Notice that this type of linear EOS is also satisfied by 
local variables: tt — w$q. The interaction is assumed 
to be proportional to the conservation term of the CDM 
dust component: 



J, = 3 e p, TL,, 



(87) 



where e is an adjustable constant parameter (to compare 
with the notation of equation (19) of [7] we have e = —a). 
For this functional dependence of J, we get 



[In J*]' _ [In p., +hiH,]' 
ZR'/R ~ 



3R'/R 



Using ( 50 1 , it is straightforward to verify that ( 87 ) and 



are consistent with J = 3 epTL in the local conserva- 
tion law ( 49 1 . The evolution equations are the following 



particular case of (51 1 
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7i* = 

P* = 
q* = 
= 

8® = 
5™ = 



-Hi- -[p* + (l + 3w )q*], 
6 

—3(1 — e) p* H*, 

-3 [(1 + w ) q* + ep*] W», 



-3W 

-HjW(l + 8M) + 



K 

6TU 



P,(SM - 6^) + (1 + 3w )q*(8M - 5^) 



(89a) 

(89b) 
(89c) 
(89d) 

(89e) 
(89f) 



which is a completely determined system depending on 
two constant free parameters: wq and e. 

A frequent motivation to consider interaction between 
CDM and DE in the numerous works in the litera- 
ture [SJ E] has been to address the so-called "coincidence 
problem" . In the FLRW solutions of [6] , the energy den- 
sities of CDM and DE are found as a function of the scale 
factor, and the authors adjust e in order to obtain a ratio 
of these energies that displays an asymptotically slowly 
varying behavior. Since we are identifying FLRW energy 
densities with the background variables p* and gv, this 
procedure involves verifying the asymptotic behavior of 
the ratio p*/q*, both expressed as functions of the scale 
factor I. Inserting p = p + q and p = woq into (64b), and 
replacing wq with wq — —\wq\, we obtain 



P* 
q* 



[p*] 



£3(1 -e) ' 

(K|-e)[g,] a -e[p,Ml-l- 3 ^°l- e >) 
(\w \ - e) PQ-M 



(90a) 



so that for \wq\ > e and f > 1 we have 



P* 

q* 



g-3(\w \-e) ^ 



(90b) 



(91) 



Hence, for e > smaller than but sufficiently close to 
I wo I, the ratio p*/g* would decay very slowly. However, 
we will follow here a different motivation. 



The conservation equation ( 49 1 in its form ( 89b I can 



also be understood as describing an interaction in which 
total CDM particle numbers increase (e > 0) or decrease 
(e < 0) HI [46] . These numbers would be conserved 
if e = 0, leading to a non-interactive mixture. If we 
consider the case with spatially flat hypersurfaces 3 T 
(see Appendix B), which follows as the particular case 
[k*]i = K = for all r in and (59), the proper spa- 
tial volume associated to these spatially flat metrics leads 
then exactly to the same integral used in the definition 
of QL variables ((22). Thus, the total number of CDM 
particles inside a comoving sphere (centered at r = 0) in 
any given hypersurface 3 T is 



N = 4tt 



t , 47T t 

iR 2 R'dx = —n*R 3 , 



(92) 



where n — p/mo is the particle number density for CDM 
particles of rest mass mo. With the help of (61 1 and 
(89b I, and bearing in mind that n* = p^/m , equation 



(92 1 leads to the following particle number scaling law: 

N = N i e 3e . (93) 

where A^(r) is an initial particle number distribution. 
Since at earlier stages the universe must have been dom- 
inated by CDM while presently DE dominates, it is rea- 
sonable to consider a scenario with CDM particle num- 



bers decaying and so ( 90 1 and ( 93 1 suggest considering 
e = — |e| < 0, but with |e| < |w[7lnstead of the "coin- 
cidence problem" we consider how the dynamics of the 
mixture is affected by decaying numbers of CDM parti- 
cles in a late time universe. 

For the numerical examination of the evolution equa- 
tions (89 1, we consider for the spatially flat case [k*]i = 
the coordinate choice (63 1 with / = tanr, and the pa- 
rameter values: w = —0.9 that is close to a cosmological 
constant, and e = —0.1, so that CDM particles decay. 
Appropriate initial value functions that comply with an 
evolution free from shell crossings (satisfying with condi- 
tion (67l) are given in Appendix D. 

The central and asymptotic values given by the initial 
conditions in (110 1 (Appendix D) are adequate for den- 
sity clumps in a matter dominated era. Hence, we can 
take hj 1 as an early universe Hubble scale factor. In 
order to verify that the evolution of the mixture yields 
CDM and DE energy densities that are minimally com- 
patible with observations, we define the following Omega 
parameters: 



~3«r 

AW EE ±* 



(94a) 
(94b) 



whose time evolution is plotted in figure [5] for the range 
< r < 7r/2. As shown by this figure, CDM is initially 
dominant with Q( p ' « 1 and 0,^ ss at the initial time 
t = 0, while as r proceeds Cl^ -> and fjW _> 1. if 
instead of (94 1, we use Omega parameters defined with 
local quantities, p, q and Ti., we obtain almost identical 
curves. From figure [5] we can roughly identify present 
cosmic era values around r = 1.25. 
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FIG. 5: The Omega parameter for CDM and DE. The 

figure displays the time evolution for the full radial range 
of the parameters fi( p ' (top panel) and Cl^ (bottom panel) 



defined by (941. We can roughly identify the present day era 
as t — 1.25. 



FIG. 6: Relative fluctuations of CDM density. The 

figure displays the relative fluctuation of CDM, <r p ' . This 
function passes from negative to positive values for all r, in- 
dicating density clumps evolving towards density voids, as 
shown in figure 9. 




In order to examine the density profiles of CDM and 



DE it is useful to recall from ( 23a ) that 



§ (p) 
6® 



1 



p*R 3 ju 
1 



q*R 3 Jo 



p' R 3 dx, 
4 R 3 dx, 



(95a) 
(95b) 



so that (for positive and q*) the sign of the relative 
fluctuations determine whether p and q present radial 
profiles of over-densities or voids in a domain of r around 
a symmetry center. In particular, it is possible to provide 
initial conditions in which an initial over-density evolves 
into a void. This situation has been reported in dust 
sources [101 [S3], and in particular, initial conditions (110) 
yield such an evolution for the CDM densities and 
The relative fluctuation d^ p ' is depicted by figure 
showing how for all layers this function passes from neg- 



FIG. 7: Relative fluctuations of DE density. The figure 
displays the function <5' 9 '. This function takes small nega- 
tive values (~ 1CP 2 ) for all layers, indicating low depth den- 
sity clumps evolving towards homogeneity. Since pressure 
anisotropy is the fluctuation of DE pressure, it is proportional 
to and so we have V/p ~ 10 -2 . 



ative to positive as r evolves. Hence, from (95a I, we have 



an initial CDM clump evolving into a configuration with 
radial profiles of void. Notice that near the center r = 
and the asymptotic range r — > n/2 we have i5' p ' close 
to zero, nevertheless, it passes from negative to positive 
values for all layers. This "clump into void" evolution is 
corroborated in the top panels of figure [8] which displays 
for various values of r the radial profiles of p and p*, 
normalized to their central values. As revealed by this 
figure, it is striking to see how an initial CDM density 
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PJ p,\ 




iJli 



FIG. 8: Radial profiles of CDM and DE energy den- 
sities. The panels display the QL and local densities for 
CDM (top two panels) and DE (bottom panels), normalized 
by their central values, as functions of r for the following con- 
stant values running from r = to r = 4.0 in intervals of 
0.2, from bottom to top. The curves show how local and QL 
variables exhibit the same qualitative behavior. The figure 
confirms that an initial CDM density clump evolves towards 
a deep density void, while DE shallow clumps slowly tend to 
homogeneity. 



clump profile evolves into a clear and deep void. The 
values of r start at t — and go up to r = in steps of 
0.2, so if we take r ss 1.25 as indicative of our present era 
(figure[5|, then the top panels of figure [8] predict that the 
CDM radial profile in our era is a very shallow clump, or 
in other words, that we would be in the middle of the 
transition from clumps into voids. 

On the other hand, as shown by figure [7] the relative 
fluctuation of DE, S^ q \ is negative and small (~ 1CP 2 ) for 
all layers and for all r. Hence, the radial DE profile is that 
of shallow clumps that slowly tend to homogeneize. This 
evolution is confirmed by plotting in the lower panels 
of figure [8] the radial profiles of DE energy density q, 
(normalized by their central values) for various values 
of constant r. From (30b I, ( |50| and (86 1 the ratio of 
anisotropic to isotropic pressures is 2V~pK — 5^/(1 + 
5^) w ~ 10~ 2 , and so this anisotropy is very small. 

While the 'clump becoming void' evolution that we 
have found (see figures [6] and top panels of figure [8]) is 
still present even when e — (no interaction, conserved 
N) , the effect of a larger (more negative) e is to increase 
the depth of the resulting voids and the speed at which 
they occur. This is consistent, because as CDM particle 
numbers decay faster one would expect the voids to be 
deeper and to arise faster. At the same time, the DE 



energy density profile becomes homogeneous at a faster 
rate. However, if we take e > 0, so that CDM particles 
are created, then we get the opposite effect: the DE pro- 
files pass from clumps to voids and CDM has the profile 
of clumps. We have used here very idealized parameters 
and do not claim that the model mixture is "realistic" , 
but by following a more careful and comprehensive ap- 
proach we could consider this type of models as describ- 
ing mechanisms to explain the relation between genera- 
tion of voids and an interaction involving the annihilation 
of CDM particles into DE. We will examine this matter 
in a separate article. 



XIII. CONCLUSION 

We have presented a formalism, based on quasi-local 
(QL) variables, describing a large class of spherically 
symmetric models (LTB spacetimes) as non-linear, gauge 
invariant and covariant perturbations of a FLRW formal 
background. The integration of the resulting evolution 
equations require numerical work, but this task can be 
handled with relatively simple numerical techniques. A 
summary of how this formalism was motivated, intro- 
duced, developed and applied is furnished in the Intro- 
duction (section II). 

Since LTB spacetimes are compatible with a wide va- 
riety of "equations of state" (EOS) and physical as- 
sumptions, these models and the formalism developed 
for them are useful theoretical tools to examine a large 
amount of cosmological sources of interest that have been 
studied only in a FLRW context (or linear perturba- 
tions) HJ|3|31H|51|B1|7IIE1|S]. While the inhomogene- 
ity of LTB spacetimes is certainly idealized, the resulting 
models are non-trivial and still exhibit non-linear behav- 
ior, and can be useful to examine important phenomena 
(gravitational collapse, void formation) that cannot be 
studied with FLRW models or their linear perturbations. 



As we argued in section VII the EOS that deter- 
mines the evolution equations is given between the QL 
energy density and pressure, so that local density and 
pressure and pressure are related by means of correla- 
tion terms that could contain non-local information that 
could be important for self-gravitating sources. Since 
the fundamental nature of dark matter and dark en- 
ergy is not known, we cannot rule the possibility of 
non-local effects playing an important role when we ex- 
amine these sources under inhomogeneous conditions, 
at least in a scale smaller than 300 Mpc. We be- 
lieve that LTB models can be useful to examine these 
and other non-linear effects, not only for dark matter 
or dark energy, but also in theoretical proposals that 
seek to explain cosmic acceleration without dark en- 
ergy [^|ISlini|ig|Tg||2a|2IJ|22||23|22|23|3g. As 
shown elsewhere [25| [27], the formalism presented here 
can be applied to understand the issue of "back-reaction" 
that appears in various averaging formalisms [THl HH1 HI] 
aiming to explore this type of alternative explanations. 
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An important result of this article is the smooth and 
fully relativistic generalization of the Newtonian "top 
hat" models (section XI I , which are useful toy models for 
structure formation and are widely used in Astrophysical 
literature. One of the numerical examples of application 
of the formalism that we provided (section XII A I is a 
Chaplygin gas "top hat" model that describes the for- 
mation of a local black hole smoothly embedded in an 
expanding Chaplygin gas FLRW universe. 

The second numerical application example is that of an 
interactive mixture of cold dark matter (CDM) and dark 
energy. By interpreting the interaction term in terms 
of particle creation, we obtained a model in which initial 
clump (or overdensity) radial profiles of CDM evolve into 
void profiles as CDM decays into dark energy. While this 



model is not "realistic" , it serves to illustrate how LTB 
spacetimes can illustrate important non-linear features 
of dark energy/matter sources that cannot be studied 
within a homogeneous context (nor with linear perturba- 
tions) . As shown in [HJ |7] , mixtures of this type can be 
studied by means of a dynamical systems approach. This 
type of approach is compatible with the evolution equa- 
tions of LTB spacetimes (autonomous equations), and so 
a generalization can be readily made of these dynamical 
system studies in a FLRW context. In fact, a dynamical 
systems study of LTB dust solutions has been achieved 
using the same QL variables [3D] . The extension of this 
work for LTB spacetimes with nonzero pressure is being 
currently undertaken. 
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Appendix A: Initial conditions. 



(341 



K* = (98a) 

JC,^ - piJ^-2H^ n \ (98b) 

so that we can always eliminate either one of n*, or 
§(H) m terms of K, and 8^. Thus, an alternative 
set of initial conditions follows by selecting 



\p*]u [}C*]i, 
(1581) and (|98| as 



[P*]i = [p«]»([M*]i)> 



(99) 



and then 8j_ as in (97 1, while [H*]i and 8. 



(H) 



follow from 



1/2 



r (k/3) - [£,]', 



r(W) 



6 (k/3) - [£. 



(100) 



We have used initial conditions given by (|99j) and (100 1 
in the numeric examples of section |XII| 



One possibility for setting up the initial value functions 



for the system ( 36 1 for an EOS ( 40 1 is simply to select 



functional forms for 



\p*]i, [7U]i, \p*]i = [p*]i([M*]i), 



(96) 



(subscript i denotes evaluation at ti — 0, or r = 0) and 
then, bearing in mind that 1^ = 1, and using ( |58[ ) to 
obtain 



r \fi 
3 [n 



I 

*Jz 



3 [H*]i 



(97) 



However, a more intuitive choice involves considering as 
an initial value function the scalar curvature [fc*]j = 
[/C*]j//ig, since there is a connection between the kine- 
matic evolution of the scale factor I (or V) and the sign 
of this curvature. This follows from the Friedman equa- 
tion (64c I, so that > is a necessary condition for 
= I — H.*, indicating a "bounce" of £(t, r). On the 
other hand, [fc*]j < is a sufficient condition for ^ being 
a monotonous function of t. 

The fact that the sign of [k*]i determines the kinematic 
behavior of I is analogous to the relation of the scale 
factor a to the constant "curvature index" ko in FLRW 
spacetimes. The difference is that in FLRW spacetimcs 
fco is a constant and thus all comoving layers in a given 
spacetime have the same kinematic evolution, while in 
LTB spacetimes [fa*]i = [k*]i(r). Therefore, if this func- 
tion changes sign in its domain of r various comoving 
layers in the same LTB spacetime can have different kine- 
matic evolution (as in the example of the Chaplygin gas 
in section 



XII A I 



In order to set initial conditions in terms of JC, /C* and 
, we note that these curvature scalars can be obta ined 

/'* 



7i*, 5^, from the Friedman equation (28 1 and 



Appendix B: Special initial conditions: the zero 
curvature case. 



It is evident from equations ( 62 1 and ( 98 ) that the 



following choice of initial conditions: 

[/CJi = 8^ =0, =► /C* = 8 {K) = 0, (101) 



leads to a simplified and important particular case that 
can be identified as LTB spacetimes with zero spatial 
curvature. Conditions (11011) lead to the constraints: 



Hi = 



5W = 



(102) 



which allows us to eliminate a nd 8 ^ in terms of /j* 
and S^' . The evolution equations ( 36 1 then simplify con- 



siderably and become the following system of two equa- 
tions: 



3/2 



= — v3k/x 

V3K 1/2 



[1 + w], (103a) 

(103b) 



where w = p*/fj,* and there is no restriction on the 
EOS. These equations can be integrated analytically for 
many choices of EOS (see [131111]), but, in general, they 
also need to be integrated numerically. Initial conditions 



( 101 1 were used in the mixture example of section XII B 



Appendix C: Initial conditions for the Chaplygin 
gas. 

Dimensionless forms for the initial value functions 
[fj.*]i, [p*]i, [Tt*]i for the Chaplygin gas "top hat" model 
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of section XII A follow from applying the Chaplygin gas 
EOS to (|58| and (|100[), and using the prescription ([74]) 



and the dimcnsionlcss parameters A and m, defined in 
(82 1. The resulting initial value functions are: 



(/<■) 



3ft 2 



K + A 2 ] 1/2 , 
A 2 



An) 



3ft, 2 : [ to; 

3 n<*° ' 

3 [w. 



+ A 2 l 1 /2- 



[fc*] 



6 n^-fc 



(104a) 
(104b) 
(104c) 
(104d) 
(104e) 



where fej = [/C*]i/ft 2 



The functions in (104 1 must lead to solutions of the 



evolution equations complying with the matching and 
regularity conditions ( 69 1— ( 72 ) and with the condition 



to avoid shell crossings ( 67 1 . Initial value functions satis- 



fying the above mentioned requirements follow by taking 

A = 0.001 and assuming that f2j and [k*]i are given by 
the simple polynomial ansatz 



Z = 



z c 



(15 - 24a; + 10x 2 ) (Z c - Z b ) , 



< x < 1 
x > 1 

(105) 

where x = r/rf,, so that Z c = Z(0) and Zb = Z(l). 
Central and background values are given by 



n^\o) = i.i, n^(i) = i.o, 
[fcJ^O) = o.i, [fe.],(i) = o, 



(106a) 
(106b) 



corresponding to an initial overdensity with small density 
contrast and with small positive curvature matched to a 
spatially flat background. 



This choice allows us to use a finite coordinate range 
0<r<7r/2to examine the full asymptotic range along 
the 3 T associated with i?j — * oo. 

Since we considered also the spatially flat case fcj = 
[JC*]j = 0, then (100 1 becomes 



[7U 



3 ([M* 



1/2 



cos r sin r (k/3) - ■) 

6 («/3) ([/**]<- [g,]i) ' 



(109) 



An evolution free from shell crossings (complying with 
( 67 1 ) is obtained with dimensionless initial value func- 



tions given by the simple ansatzes 



3ft 2 
3ft 2 



= 1.2 



= 0.1 



0.6 



+- 9 tan 4 r ' 
0.05 



1 



tan 2 r 



(110a) 
(110b) 



so that the central and asymptotic values for «;[/i >t ] i /(3/i 2 ) 
are, respectively 1.8 and 1.2, while the corresponding val- 
ues for Av[<7*]i/(3/i 2 ) are 0.15 and 0.1. Also, we considered 
the value w = —0.9 that is close to a cosmological con- 
stant, and e = —0.1, so that CDM particle numbers de- 
cay. 

The remaining initial value functions, [7^] ^ S^ p \ 6^ 



Chi \ 

and 8\ , follow from relations like ( 97 1(1 00[) , but con- 
sidering that now [/C*]j = and = [p7\ 
also the coordinate choice (1071, 1108a I and (|108b|. Ini 



and 



tial conditions ( 110[ ) yield an evolution free from "shell 
crossings" (condition (67 1 is fulfilled). 



Appendix D: Initial conditions for the interactive 
CDM-DE mixture. 



In the example of the interactive DM-DE mixture of 
section |XIIB we chose initial value functions for ( 89 1 by 



considering (instead of ( 56 )) the following choice of radial 
coordinate: 



Ri = h 1 tan r, 



so that (|57b| and (|58| become 

= 1 + cos r sin r 



r = 1 



6W = 



e/i 



R'JRi 

1 A^/A, _ cos t sinr A' t 
3T R'JR, = 3T A 



(107) 

(108a) 
(108b) 



